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. Abstract 

The products of the electron width of the J/ift meson and the branching fraction of its decays to the lepton pairs were measured 
using data from the KEDR experiment at the VEPP-4M electron-positron collider. The results are 
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■■ 0.3323 ± 0.0064 (stat.) ± 0.0048 (syst.) keV, 
0.3318 ± 0.0052 (stat.) + 0.0063 (syst.) keV. 
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0.6641 ± 0.0082 (stat.) ± 0.0100 (syst.) keV, 
1.002 ± 0.021 (stat.) ± 0.013 (syst.) 
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can be used to improve the accuracy of the leptonic and full widths and test leptonic universality. 

Assuming ey. universality and using the world average value of the lepton branching fraction, we also determine the leptonic 
Q\ T H = 5.59 + 0.12keV and total T = 94.1 + 2.7keV widths of the J/ip meson. 



Keywords: J I if/ meson, leptonic width, full width, leptonic universality 
PACS: 13.20.Gd, 13.66.De, 14.40.Gx 
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Cu 1. Introduction 

_ _ i 

The // uj meson is frequently referred to as a hydrogen atom 
for QCD. The electron widths T ee of charmonium states are 
rather well predicted by potential models U Q]. The accu- 
racy in the QCD lattice calculations of T ee gradually approaches 
the experimental errors ^ . The total and leptonic widths of a 
hadronic resonance, E and Fa, describe fundamental properties 
of the strong potential . 

In this paper we report a measurement of the product of the 
electron width and the branching fraction to an e + e~ pair for 
the J/ifr meson, T ee x r ee /T. An experimental determination 



of F ee x F f ,/ 1 requires scanning the beam energy and measur- 
ing the cross section. In contrast to a measurement of the lep- 
tonic width itself, in this case knowledge of the efficiency for 
hadronic decays does not contribute to the final uncertainty. The 
problem considered can be reduced to measuring the area under 
the resonance curve for the process e + e~ — * J/ifr — > e + e~ . Ad- 
ditionally we have measured the product of the electron width 
of the J /if/ meson and the probability of its decay to the 
pair, T ee x Tfju/F. Given independent data on the branching 
fraction S ee J5[], we use this result to evaluate the leptonic Ta 
and total Y widths. 
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Figure 1 : VEPP-4M/KEDR complex with the resonant depolarization and the 
infrared light Compton backscattering facilities. 



2. VEPP-4M collider and KEDR detector 

The VEPP-4M collider [6] can operate in the broad range 
of beam energies from 1 to 6 GeV. The peak luminosity in the 
J/ift energy range is about 2 x 10 3l) cm _2 s _1 . 

One of the main features of the VEPP-4M is a possibility 
of precise energy determination. The resonant depolarization 
method 

fill 

was implemented at VEPP-4 at the very begin- 
ning of experiments in early eighties for the measurements of 
the J /if/ and ip(2S ) mass with the OLYA [9] detector and T fam- 
ily mass with the MD-1 |9[] detector. 

At VEPP-4M the accuracy of energy calibration with the 
resonant depolarization is improved to about 10~ 6 . Between 
calibrations the energy interpolation in the J/i// energy range 
has the accuracy of 6 ■ 10" 6 (^lOkeV) d. 

To monitor beam energy during data taking the infrared 
light Compton backscattering is employed (with 50-^70 keV 
precision in the J/ijj region), which was first developed at the 
BESSY-I and BESSY-II synchrotron radiation sources ifTH [rj . 

The KEDR detector 1 1311 includes the vertex detector, the 
drift chamber, the scintillation time-of-flight counters, the aero- 
gel Cherenkov counters, the barrel liquid krypton calorimeter, 
the endcap Csl calorimeter, and the muon system built in the 
yoke of a superconducting coil generating a field of 0.65 T. The 
detector also includes a tagging system to detect scattered elec- 
trons and study two-photon processes. The on-line luminosity 
is measured by two independent single bremsstrahlung moni- 
tors. 



3. Experiment description 

A data sample used for this analysis comprises 230 nb -1 
collected at 1 1 energy points in the J /if/ energy range. This cor- 
responds to approximately 15000 J/iff — > e + e~ decays. During 
this scan, 26 calibrations of the beam energy have been done 
using resonant depolarization. 

The primary trigger signal was provided by a coincidence of 
two non-adjacent scintillation counters or an energy deposition 
in the endcap calorimeter of at least 100 MeV. A veto from the 
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Figure 2: Observed cross section of e + e — » hadrons in the J/ip scan. 

endcap-calorimeter crystals closest to the beam line was used 
to suppress the machine background. 

The secondary trigger required at least two tracks in the drift 
chamber or at least one track and an energy deposition in the 
calorimeter of at least 70 MeV and the coincidence of two non- 
adjacent scintillation counters. 

The hardware triggers use the analogous output of the calorime- 
ter with reduced energy resolution. During the offline analy- 
sis real and simulated events pass through the software event 
filter which recalculates the trigger decision using a digitized 
response of the detector subsystems. The calorimeter energy 
thresholds in the event filter are toughened by a factor of 1.5 
with respect to the instrumental values, suppressing the uncer- 
tainty in the latter and their possible instability. 

Single bremsstrahlung and e + e~ — > e + e~ events at polar an- 
gles in the range between 18° and 31° (the endcap calorime- 
ter) were used in the relative measurement of luminosity. In 
order to evaluate T ee x r ee /F, it was unnecessary to measure 
the absolute luminosity. Since e + <?~ — > <? + e~ events analyzed 
here include both events from the resonance and a well known 
non-resonant QED background, it is possible to perform an ab- 
solute calibration of the luminosity along with the derivation of 

Egg X Y ee j E. 

Figure|2]shows the observed cross sections of e + e~ — > hadrons 
in the energy range. These data were used to fix the reso- 
nance peak position and to determine the beam energy spread. 
The value of the J/tfr mass agrees with the earlier VEPP-4M/KEDR 
experiments lioll . The accuracy of the energy spread was about 
2%, including variations associated with the beam current. 

4. Theoretical e + e~ -* t + (~ cross section 

The analytical expressions for the cross section of the pro- 
cess e + <?~ — > ( + €~ with radiative corrections taken into ac- 
count in the soft photon approximation were first derived by 
Ya. A. Azimov et al. in 1975 [ 14]. With some up-today modifi- 
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cations one obtains in the vicinity of a narrow resonance 
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Here W is the center-of-mass energy and Flo represents the vac- 
uum polarization operator with the resonance contribution ex- 
cluded. The terms proportional to Im T and Re T describe the 
contribution of the resonance and the interference effect, re- 
spectively. The definition of leptonic width in Eq. (Q]i- © im- 
plicitly includes vacuum polarization as recommended by PDG: 
Tex = r1L/|l - riol 2 , where Y° ee is the lowest-order QED value. 

The function T in eq. (f3]l appears from the integration (see 
lfl6ll . where one can also find the definition of T for the rela- 
tivistic Breit-Wigner amplitude) and differs from that in Ref. (l~ 
by the Trp/sinnf) factor. 

For the e + e~ final state one has 
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The goal of this analysis is a measurement of Y ee x Y ee j Y 
and F ee x T^l Y contained in the resonant terms. The precision 
of these terms in formulae ([TJ and (0 is better than 0.2%. It was 
estimated by numerical calculations beyond the soft photon ap- 
proximation according to Ref. II 1 511 . Although the interference 
terms could allow a direct measurement of Y ee (e + e~ — > e + e~) 
and JT^eY^t (e + e~ — > ii + [r), in our case we are limited by the 
statistical accuracy and theoretical uncertainty. 

To compare experimental data with the theoretical cross sec- 
tions ([1) and (0, it is necessary to perform their convolution 
with a distribution of the total beam energy which is assumed 
to be Gaussian with an energy spread crw- 
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■ exp 
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where Wo is an average cm. collision energy. 



Since the energy spread crw — 0.7 MeV is much larger than 
the intrinsic width of the J /if/ meson, the uncertainty of the cross 
section due to the knowledge of the latter is suppressed. We use 
the value Y * 0.093 MeV Jj . 

For simulating the nonresonant contribution ctq ED we use 
the calculations of II17L as implemented in two independent 
generators BHWIDElEirand MCGPJ Q . 

The resonant and interference cross sections were simulated 
using simple generators with proper angular distributions. In 
this case the initial state radiative corrections are already taken 
into account by the expressions ([TJ and (0. These formulae 
implicitly involve the branching ratios Yu /Y = S(C( n y) with the 
arbitrary number of soft photons emitted. Actual event selec- 
tion criteria can not be 100% efficient for events with additional 
photons, therefore the final state radiation must be simulated 
explicitly. This was done using the PHOTOS package [21]. 

5. Data analysis 

In our analysis we employed the simplest selection criteria 
that ensured sufficient suppression of multihadron events and 
the cosmic-ray background. The following requirements were 
imposed for e + e~ — > e + e~ events selection: 

1 . An event should have exactly two oppositely charged tracks, 
each originating from the beam intersection region, hav- 
ing a continuation in the calorimeter, and lying in the 
range of angles between the particle and beam axis from 
30° to 150°. 

2. The energy deposited in the calorimeter for each particle 
should be higher than 0.7 GeV, and the sum of the ener- 
gies of the two particles should be higher than 2 GeV. 

3. The energy deposited in the calorimeter and not associ- 
ated with the two particles considered should not exceed 
5% of the total energy deposition. 

4. The angle between selected particles should be larger than 
140° and acoplanarity less than 40°. 

Requirements for selecting e + e~ — > events are: 

1. The same as tracking criteria for e + e~ — > e + e~ . 

2. The energy deposited in the calorimeter for each particle 
should be higher than 60 MeV and less than 500 MeV, 
and the sum of the energies of the two particles should 
not be higher than 750 MeV. 

3. The energy deposited in the calorimeter and not associ- 
ated with the two particles considered should not exceed 
30% of the total energy deposition. 

4. The angle between selected particles should be larger than 
170° and acoplanarity less than 15°. 

5. The momentum for each particle should be higher than 
500MeV/c, and the sum of the momenta of the two par- 
ticles should be higher than 2 GeV/c. 

6. There is at least one time measurement in the time-of- 
flight system. The uncorrected measured time should be 
within the [-3.75 -f 10.] ns range from the beam intersec- 
tion time. 
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Figure 4: Fits to experimental data for the process e + e — > e + e in the J /if/ 
energy range for four angular ranges. 



Figure 3: Cross section of the process e + e~ — > e + e~ as a function of the electron 
scattering angle at the J/if/ peak. The points represent experimental data. The 
histograms correspond to a simulation: the dashed line represents the contribu- 
tion of the J/ift resonance,the dashed and dotted line represents the contribution 
of Bhabha scattering and the solid-line histogram is the sum of the first two. 



Figure [3] shows the distribution of selected e + e~ — > e + e~ 
events with respect to the electron scattering angle. The dis- 
played points represent the experimental values, while the his- 
tograms correspond to the simulation. At small angles Bhabha 
scattering prevails, while at large angles events of resonance de- 
cay are dominant. The interference effect is not shown since the 
presented data correspond to the J /iff peak, where the interfer- 
ence vanishes. 

In order to measure the resonance parameters, the set of 
events was divided into ten equal angular intervals from 40° to 
140°. At the z'-th energy point E, and the j-th angular interval 9j, 
the expected number of e + e~ — > e + e~ events was parameterized 

as 



N exp (Ei,8j) =S £ xX(£i)X 



cr*r r (£^./)-er(£^ ; )+ 



inter v " J> n 



r (Ei,ej)+ 



(6) 



<abha(^.^')-<abha(^.^)). 



where -£(£;) is the integrated luminosity measured by the lumi- 
nosity monitor at the z'-th energy point; cr*™ 1 ', cr*™ r and crSS, 
are the theoretical cross sections for resonance, interference and 
Bhabha contributions, respectively, sf™, Sj™ r and Sg™ bha are 
detector efficiencies obtained from simulation. The efficiencies 
differ mainly due to difference in radiative corrections. Unlike 
the Bhabha process, the initial state radiation for the narrow res- 
onance production is strongly suppressed, thus the events are 
more collinear. 

In this formula the following free parameters were used: 

1 . the product F ef x T ee l T, which determines the magnitude 
of the resonance signal; 

2. the electron width F ee , which specifies the amplitude of 
the interference wave; 

3. the coefficient %£, which provides the absolute calibra- 
tion of the luminosity monitor. 



We note that the coefficient %£ partially takes into account a 
possible difference between the actual detection efficiency and 
simulation in the case where this difference does not depend 
on the scattering angle or the beam energy (or the data taking 
time), thus a substantial cancellation of errors occurs. 

The r ee value obtained from the fit to the data has large 
statistical and systematic uncertainties caused by the smallness 
of the interference effect and the low accuracy of theoretical 
evaluation. 

Figure|4]shows our fits to the data for four angular bins. For 
this fit^ 2 /ndf = 53.7/41 taking into account only the statis- 
tical errors and^ 2 /ndf ^ 40.5/41 after converting the energy 
determination uncertainty to the cross section error. 

The joint fit in ten equal bins from 40° to 140° produces the 
following basic result: 

r ee x T e J r = 0.3323 + 0.0064 (stat.)keV, 

n £ = 93.4 ±0.7 (stat.) %, (7) 

r ee =5.7 ±0.6 (stat.) keV. 

Due to different angular distributions for Bhabha scattering 
and resonance events, subdivision of the data into several angu- 
lar bins reduces the statistical error for T ee xF ee / F by 40H-50 %. 
Here T ee has a statistical error of about 10% and agrees with the 
world average value. The same value can be obtained with a 
much higher precision using Ta X T^/r and an independent 
measurement of the branching ratio J/\j/ — > 

Similarly to ©, the expected number of e + e~ — > fi + fi~ 
events was parameterized in the form: 



ATexp(£,) =7?x x £(£,)* 

(o^*(E t ) ■ eC(Ed+ 



theoi /c-A . sim 
u '""v \^'> 6 nte 



(8) 



inter 
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< eOT (£/) • e\ 



+ F C . 



with the same meaning of %£ and £,(Ei) as in (|6). £(Ei) is mul- 
tiplied by the sum of the products of theoretical cross sections 
for resonance, interference and QED background and detection 
efficiencies as obtained from simulated data, fi^ was fixed from 
the result (0 and T, is the live data taking time. Unlike ©, there 
is only one angular bin from 40° to 140°. 
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Figure 5: Fit to experimental data for e + e~ — > /J + /i process in the J/if/ energy 
range. 

The following free parameters were used: 

1 . the product Y ee x Y Mfl / Y, which determines the magnitude 
of the resonance signal; 

2. the square root of electron and muon widths ^JT\e^^, 
which specifies the amplitude of the interference wave; 

3. the rate of cosmic events, f cosm i c , that passed the selec- 
tion criteria for the <? + e~ — > y?pT events. 

Due to variations of luminosity during the experiment it is pos- 
sible to separate the contribution of cosmic events (f cosmic ■ Ti) 
from that of the nonresonant background ((X? eor (.Ej) ■ s^™(Ei) ■ 

X( /•:,)). 

Figure [5] shows our fit to the e + e~ — > [T data. It yields 
the following result: 

Y ee xyr = 0.3318 ± 0.0052 (stat.) keV, 

/ (9) 

JYee x r w = 5.6 ± 0.7 (stat.)keV. 

As can be seen from ((9), the statistical error of T ee X Y^j Y is 
about 1.6%. 

6. Discussion of systematic uncertainties 

The most significant systematic uncertainties in the Y ee X 
r<, e /r and Y ee X T^/T measurements are listed in Tables [1] 
and [2] respectively. A few dominant sources of uncertainty are 
briefly described below. 

A rather large uncertainty of 0.8% common for the electron 
and muon channels is due to the luminosity monitor instability. 
It was estimated from comparing the results obtained using the 
on-line luminosity of the single bremsstrahlung monitor and the 
off-line luminosity measured by the e + e~ scattering in the end- 
cap calorimeter. 

The essential source of uncertainty is an imperfection of the 
detector response simulation resulting in the errors in the trigger 
and offline event selection efficiencies. 

To correct the offline event selection efficiency, two high- 
purity control samples of e + e~ events were prepared. The first 



Table 1: Systematic uncertainties in T ee X r ee /r. 



Systematic uncertainty source 


Error, % 


Luminosity monitor instability 


0.8 


Offline event selection 


0.7 


Trigger efficiency 


0.5 


Energy spread accuracy 


0.2 


Beam energy measurement (10-30keV) 


0.3 


Fiducial volume cut 


0.2 


Calculation of radiative corrections 


0.2 


Cross section for Bhabha (MC generators) 


n a 
U.4 


rinai state lauiation (rriuivJo) 


n a 


Background from J/ift decays 


U.z 


Fitting procedure 


n i 
U.z 


lotal 


1 A 


Table 2: Systematic uncertainties in Y ee X T^j 


r. 


Systematic uncertainty source 


Error, % 


Luminosity monitor instability 


0.8 


Absolute luminosity calibration by e + e~ data 


1.2 


Trigger efficiency 


0.5 


Energy spread accuracy 


0.4 


Beam energy measurement (10-30keV) 


0.5 


Fiducial volume cut 


0.2 


Calculation of radiative corrections 


0.2 


Final state radiation (PHOTOS) 


0.5 


Nonresonant background 


0.1 


Background from J/tfr decays 


0.6 


Total 


1.9 



sample selected using the LRr-calorimeter data only was em- 
ployed to determine the tracking system efficiency, the second 
sample obtained using mostly the tracking system data allows 
one to check calorimeter related cuts. Each sample contains 
about 70% of all events used in the analysis. The same analy- 
sis was performed with simulated data. The corrections already 
taken into account in (|7]) were 

6 Y ee xYee/Y = 0.8 ± 0.6 (stat.) ±0.4(syst.)%, 
6R £ = 1.7 +0.5 (stat.) + 0.5 (syst.) %. 

The statistical error of the efficiency determination is approx- 
imately three times less than that of the final result due to the 
binomial distribution in the number of lost events. The residual 
systematic error is due to an incomplete event sample employed 
for the correction and the efficiency difference for the resonance 
decays and the continuum events. The variation of %£ is greater 
than that of the main result illustrating the cancellation of un- 
certainties mentioned in Sec. [5] 

Three contributions dominate the trigger efficiency uncer- 
tainty. The inefficiency of the time-of-flight counters used in 
the first level trigger was studied using the cosmic ray events 
and equals 0.3%. The second contribution comes from the cut 
on the number of the vertex detector tubes hit in the event. It 
was used in the software trigger level for the machine back- 
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ground suppression. Some fraction of events was accepted un- 
conditionally to check the cut. The third contribution is due to 
the veto from the Csl crystals nearest to the beam line. It is 
negligible for the resonance decays and reaches 0.4% for con- 
tinuum events for which the initial state radiation is not sup- 
pressed. The quoted value was obtained varying the threshold 
in the event filter within its uncertainty. The accidental signal- 
background coincidences were taken into account by the veto 
rate with much better accuracy. 

The uncertainty of the theoretical Bhabha cross section was 
estimated comparing the results obtained with the BHWIDE [ 19] 
and MCGPJ B20I1 event generators. It agrees with the accuracies 
of the generators quoted by the authors. 

The dominant uncertainty of the Y ee X Y^/ T result is as- 
sociated with the absolute luminosity calibration done in the 
e + e~-channel. It includes the accuracy of the Bhabha event gen- 
erators, the statistical error of %£ from © and the residual ef- 
ficiency difference for e + <?~ and ^ yT events after a correction 
using simulated data. 

To determine this residual difference and make a proper cor- 
rection, samples of real and simulated quasi-collinear events 
were selected using an alternative track reconstruction code find- 
ing a single track with a kink at the point of the closest approach 
to the beam line. Then the standard analysis procedure was per- 
formed for these samples yielding the double ratio 



„ es P\ 



/e exp \ 



1.005 + 0.005 (stat.) + 0.008 (syst). (11) 



The sample selected contains about 80% and 50% of all ytyT 
and <? + e~ events, respectively. The systematic error of the ratio 
reflects the incompleteness of the samples. 

The trigger veto uncertainty is the same as for J /if/ — > e + e~ 
decay. 

The background for J/iJ/ — > /i 4 yT decay from hadronic de- 
cays of J/tf/ was estimated with the help of the muon system. It 
contributes 1 .5 ± 0.6 % to the selected yfyT events. The estima- 
tion agrees with the simulation results. This correction as well 
as the correction ( fTTT ) have already been taken into account in 
the r ee xr w /r result©. 

Due to the high precision in the energy determination by 
the resonant depolarization method 12211 . the corresponding er- 
rors (peak position, energy spread, and energy at a point) are 
relatively small. 

In calculating the cross section for resonance production 
with formulae © and © we used the PDG value of the total 
width F from [5J. Its error is about 2%, which gives a ~ 0.05% 
contribution to the error in our result. 

The fiducial volume cut 40° < 6 < 140° was applied using 
the tracking system and the strip system of the LKr calorimeter. 
The difference of results provides a conservative uncertainty es- 
timate. 

All other uncertainties are rather clear. More detail can be 
found in j23ll . 

All the uncertainties for Y ee x T ee IY added in quadrature 
yield a systematic error of 1.4%. All uncertainties for Y ee x 
Tun/ T added in quadrature yield a systematic error of 1.9%. 
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Figure 6: Comparison of T ee X T ee j T and T ee X Tun/ T measured in different 
experiments mentioned in |5] with KEDR 2009 results. The grey strip is for the 
world average T ee X T MU / T value. 



7. Results and Conclusion 

The new measurement of the Y ee xY ee /Y and Y ee x Y^/ Y 
has been performed at the VEPP-4M collider using the KEDR 
detector. The following results have been obtained: 

Y ee x T ee / r = 0.3323 ± 0.0064 (stat.) ± 0.0048 (syst.) keV, 
r ee x T m l T = 0.3318 ± 0.0052 (stat.) ± 0.0063 (syst.) keV. 

Previously, Y ee x Y ee j Y was measured in the DASP exper- 
iment in 1979 [24] with a precision of about 6%. The result 
obtained in the present study improves the accuracy by a factor 
greater than two. The most precise previous measurements of 
Y ee x r^/ T were made in the BaBar [25] and CLEO-c [26] 
experiments, both with the ISR technique. 

Figure [6] shows the comparison of our results with those of 
the previous experiments. The grey line shows the PDG average 
and the error for the Y ee x Y m l Y product measurement. The 
new KEDR results are the most precise. Results are in good 
agreement with each other and with the world average value of 

Y ee x r^, / r. 

From the direct measurements of the products above one 
can extract the leptonic and full width of the resonance as well 
as test leptonic universality. For the former one should calculate 
the sum of Y ee xY ee /Y and Y ee X Y^/ Y, while for the latter the 
ratio of these quantities can be used. 

While estimating uncertainties of Y ee x (Y ee + Y^) Y and 
r ce / Y Mfl correlations between Y ee xY ee /Y and Y ee x Y Mfl / Y sys- 
tematic errors were taken into account: 

Y ee x (Y ee + T m )l Y = 0.6641 ± 0.0082 (stat.) ± 0.0100 (syst.)keV, 
T ee l Y m = 1.002 ± 0.021 (stat.) ± 0.013 (syst.) 
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In contrast to the F ee x Y ee j Y and T ee x T m f F values, the ratio 
Tec/T^u is not sensitive to the absolute luminosity calibration. 
Therefore, the parameter has been fixed in the fit and the 
relative statistical uncertainty of the Y ee j F m value is less than 

that ofr ee x(r ee + r wt )/r. 

With the assumption of leptonic universality and using in- 
dependent data on the branching fraction S(J/iff — > e + e~) = 
(5.94 ± 0.06)% S, the leptonic and total widths of the J /iff 
meson were determined: 

T ce = 5.59 + 0.12 keV, 
T = 94.1 ±2.7 keV. 

These results are in good agreement with the world average |5] 



and with the results from the BaBar 02511 and CLEO-c M26I1 ex- 
periments. 
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